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Abstract. We study the Langevin dynamics of a dipole diffusing in a random electrical field E 
derived from a quenched Gaussian potential. We show that in a suitable adiabatic limit (where the 
dynamics of the dipole moment is much faster than the dynamics of its position), one can reduce 
the coupled stochastic equations to an effective Langevin equation for a particle diffusing in an 
effective potential with a spatially varying and anisotropic local diffusivity Kij. Analytic results, 
close to the adiabatic limit, for the diffusion constant Kg, are found in one dimension and a finite 
temperature dynamical transition is found. The system is also studied numerically. In particular, 
we study the anomalous diffusion exponent in the low temperature regime. Our findings strongly 
support the conclusion that the location of the dynamical transition and the anomalous diffusion 
exponents are determined by purely static considerations, i.e. they are independent of the relative 
values of the diffusion constants of the particle position and its dipole moment. 
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1. Introduction 



Computing bulk transport properties in random media is an important physical problem ||^, 
^, ^, 1^, 0, P having many applications. Two important models of random media have been 
widely studied (i) where the randomness is due to a random potential and (ii) where the local 
diffusions constant is random. There are mathematical relations between the problems (i) and 
(ii) and problems of type (ii) are also related to a wide variety of physical problems such as the 
effective conductivity and dielectric constants of random conductors and dielectrics and the effective 
permeability of porous media [Q . An important quantity to understand those properties is the late 
time diffusion constant of a Brownian tracer particle in interaction with the medium. In the system 
we will study, the local transport properties are modified by the interaction of the dipole moment 
of the tracer particle with a random electric field drawn from an appropriate statistical ensemble. 
If the electric field disorder is statistically invariant under translation and short range correlated in 
space, we expect that the diffusion constant and mobility due to an applied force are self-averaging 
in the regime where the diffusion is normal. If the transport is normal, then in the long time limit 
{t oo), the mean squared displacement of the particle position behaves as 

((x(t) - x(0))2) - 2i?M, (1) 

where D is the dimension of space and Kg is the late time diffusion or effective diffusion 
constant. In general there are subdominant corrections to the above, which depend strongly on the 
dimensionality of the problem ||, |l^ . 

If the random field is frozen or evolves over very long time scales with respect to those of the tracer 
particle, then one is in the situation of quenched disorder, where the potential does not depend on 
time. The motivation for studying models of particles diffusing in a random potential comes from 
the fact that such systems arise very naturally in nature; for instance in zeolites where random 
electric fields are generated by the presence of frozen charged impurities [ pT| . If the particle has a 
net charge, then the particle interacts with the electric field and the case of diffusion in a quenched 
scalar field has been extensively studied in the literature ^ ^, ^ |l2|. However, if the particle 
is polarizable but without charge, it also interacts with the field. In this case, the problem is quite 



different and much less well studied |I2, We also note that systems with quenched disorder, 
spin glasses for example, are often good paradigms for systems having structural glass transition 
(where no random field is present). It is often argued heuristically, that for sufhciently complex 
and frustrated systems, a single particle in the system sees an effectively random potential due 
to the other particles. At a mean field level, there exist models, where this analogy has been 
used successfully to analyze the statistical mechanics of frustrated but non-disordered models. In 
practice one can have two models, one with quenched disorder and the other without but highly 
frustrated, which exhibit the same thermodynamics in the high temperature phase and the same 
glass transition at low temperatures ||lj, |T^, . Even if the frustrated non-random system possesses 
a crystalline ground state, not shared by the disordered system, this fact is practically irrelevant as 
this state is dynamically never attained. Therefore, these models are often adopted as toy models 
for structural glass transition. In the liquid phase, we expect the tracer particle to have a non-zero 
diffusion constant k;. Now consider the situation in which the same particle diffuses in a quenched 
background, where all the other particles have been frozen in a particular configuration. A realistic 
choice would be to select a configuration from a Gibbs-Boltzmann equilibrium ensemble. If Ke is 
the diffusion constant, it has been shown |0, uM that < ki which makes physical sense, as if the 
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background particles can move about, the cages, which trap the tracer , wiU break up on some time 
scale and free it to disperse more quickly than in the quenched case. Moreover, in some special 
cases, Ke can vanish at a critical temperature or disorder strength 12 1. If the system has a 



finite correlation length Iq (which will be the case, as our field will be short range correlated), the 
diffusion constant can be used to give an effective relaxation time r 

r^^. (2) 

If Ke — > 0, it means, that we have a diverging timescale in the same way in which a structural 
glass has a diverging timescale which, experimentally, can be extracted from the divergence of the 
liquid's viscosity. 

In this paper we study the physical case of dipoles diffusing in a random electric field E(x) which 
is spatially varying though time independent (or quenched). It is generated by a random potential 
(/)(x) which gives E(x) = — V0(x), with correlation function 

(0(x)0(x'))=A(|x-x'|). (3) 

The potential (/)(x) is thus statistically isotropic and invariant by translation in space. The most 
convenient choice is to take to be Gaussian. If the dipole moment, denoted by p, is modelled 
as two opposite charges connected to a Harmonic spring and x denotes the position of the dipole 
centre, then the total energy of a particle at the point (p, x) is given by 

i7(x,p) = i-p2-p.E(x). (4) 

The first term is the Harmonic energy of the spring, and x is simply the dipole polarisability. The 
second term is the energy of the dipole's interaction with the field E. The partition function for 
the system, where x is confined to a volume denoted by V of a Z? dimensional space, is given, up 
to an overall factor, by 

Z = [ dii [ dpexpf-/3^+/3p-E(x)y (5) 



V 



2x 



where (3 — l/kgT is the inverse temperature. If we trace over the dependence on p, we find an 
effective partition function for the variable x given, again up to a constant, by 

Ze//- / rfxexp(~/3y(x)), (6) 
Jv 

where V is the effective potential for x and is given by 

VM--'^^ (T) 

This represents the case where the dipole moment adapt instantaneously to the external field and 
motivated the study in |l^ of diffusion in non-Gaussian potentials. In [|l^ the diffusion constant 
for a particle diffusing in non Gaussian potentials such as given by equation (|^) was computed 
exactly in one dimension, and in a self renormalization group scheme was developed to study 
the problem in higher dimensions. 

In jl^, |l^ a critical temperature at which the diffusion constant vanishes, signalling a dynamical 
transition from a normal diffusive regime to a subdiffusive regime was identified. We thus expect 
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that, in a suitable adiabatic limit where the dynamics of the dipole degree of freedom p is much 
more rapid than that of the spatial variable x, that the effective potential seen by the variable x is 
V as defined by equation (Q). 

If we take an overdamped Langevin dynamics for both the dipole and positional degrees of freedom, 
the equations of motion are given by 



dxi dEj r- — 

— = "'^'■^Pl~Q^ + V^fixVxi (8) 

"'^^ - - Kpl3 f ^ - Ei) + y/2^r,p, . (9) 



dt ■^'^-Vx 

In the above equations Kx is the bare diffusion constant for the spatial variable and the other 
diffusion constant Hp sets the time scale for the relaxation of the dipole and the adiabatic limit 
is where Kp — > oo. In the dynamical equations (||) and (^, terms of the form Ai represent the 
component of the vector A in the direction i ; we have also used the Einstein summation convention 
and will stick with this convention throughout the paper. The noise terms are white noise and their 
correlation functions are given by 

iVp, {i)Vp, (s)) = {Vx, it)Vx, (s)) = 6^j6(t - s) and {rjp^ {t)Tjx, (s)) = 0. (10) 
If we rewrite (^) and substitute it into equation (^) we find 



dxi 
It 




(U) 



If we now take the adiabatic limit [np ^ Kx), equation (|Tl|) reduces to the Langevin equation for a 
particle in the potential ^(x) defined by equation (|7|). In appendix A, we will rederive this result in 
a more rigorous way and also see, at first order, the effect of a finite value of Kp. We thus expect, that 
for very large values of Kp, we find the same dynamics for x studied in [|lO|, Consequently 
there should also be a dynamical transition in this problem when the field cj) is Gaussian. We can 
now ask the question, what happens if the time scale for the relaxation of the dipole moment is 
non zero? Will the dynamical transition remain or is it a pathology of the limit Kp — s- cx)? 
In the appendix A we derive the effective dynamics for the marginal distribution of x to order Kx/ Hp- 
To this order the effective dynamics of x still can be described by a Langevin equation, with the 
same potential as equation but with a spatially varying non- isotropic diffusion constant. This 
is a rather remarkable fact and can be shown using operator projection techniques [l9[| . However 
we will present a derivation based on a direct manipulation of the Langevin equations similar to 
that of We chose this route as it gives a physical feeling for why the effective process for x 
is to this order Markovian and also because the computation in spaces of dimension greater than 
one are more straightforward within this formalism. We also compute the first order corrections to 
the effective diffusion constant in the high temperature limit via a Kubo formula for the effective 
diffusion constant. 

In section (3) we will present exact result for the diffusion constant in one dimension. In section 
(4) we will confront our results with numerical simulation (stochastic second order Runge-Kutta) 
of the coupled Langevin equations (H) and (|^) in the diffusive and sub-diffusive regime. Finally in 
section 5 we will conclude and discuss our results. 
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2. Large Kp and small /3 approximations 



The time scale for the relaxation of the dipoles should be proportional to and thus, the adiabatic 
limit (where the dipoles adapt very quickly to the local field) will correspond to the limit where Kp 
becomes large. Interestingly, in this limit, one can simplify the coupled equations (||) and (^ to an 
effective Langevin equation for x up to the order 0{kx/ Hp) by direct manipulation of the Langevin 
equations The exact derivation is a little long and technical and for clarity's sake it is thus 

given in appendix A. As a result, we can write down the following Fokker-Planck equation which 
describe the effective process for x 



— = H = — 
dt dxk 



dxi dxj 



(12) 



Where V is the effective potential given by equation (^ and is a spatially varying and anisotropic 
diffusivity tensor given by 



-X 



, dEk dEk 



(13) 



The Fokker-Planck equation (|T^ has the correct Gibbs Boltzmann equilibrium distribution with 
the effective potential V and from this one can write down a Langevin equation which corresponds 
to the process 



dxi _ ^ d xE. 
dt dxi 2 



(14) 



dxr' 

We note that if we take the limit where Kp — > oo in ([1^), the diffusivity reduce to Kij 
and we recover the Langevin equation for a particle in the potential T^(x) rigorously. In equation 
(A. 21) of appendix A we performed an expansion in a^^ (with a — KpP/x) and assumed a to 
be large. However if Kp is large but finite we see that as [3 becomes small (the high temperature 
limit) the expansion will fail. To predict the behavior of Ke at high temperature we will thus use 
a Kubo formula for the effective diffusion constant. Integrating the stochastic differential equation 
(||) between and t we have 

x^{t) - Xi{0) = V2:^Bi{t) + Pk^ f dspj{s).W^EJ{x{s)) (15) 

Where Bt is a standard D-dimensional Brownian motion with (Bj ) = 2Dt and Xi{0) is the initial 
position in the direction i. Thus squaring the above equation and taking the average yields 



2K,t = {{x,{t) - x,{0)y) + 2K,P{{x,{t) - a;,(0)) . / dspj(s).V,£;,(x(s))) 



H^^(3f / ds / ds'{p,{s).y,E,{■K{s))pk{s').y^Ek{As'))■ 
Jo Jo 



(16) 



Using the property of detail balance, which ensures time translation invariance at equilibrium, we 
have for any two functions A and B that 



(A(x(0)i?(x(.))) = (A(x(t-.))i?(x(0))), 



(17) 
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and the Onsager symmetry relation 

{A{^{t))B{^{s))) = (i?(x(t))yl(x(s))). (18) 

Assuming that the system starts in equihbrium ||2l| one can thus apply ( p^ ) and ( |l8| ) to the second 
term of ( p^ and show that it should vanish 

{{x^{t) - x,{0)) . [ dspj{s).V,Ej{x{s))) 
Jo 

t 

ds[{x^{t)pj{s).W.,Ej{x{s))) - {x,{0)pj{s).W^Ej{^{s)))] 



t 



= / ds [{x,it - s)p,(0).V,£;,(x(0))) - (a:,(0)pj(s).V,£;,(x(s)))] 
Jo 

= f ds[(p,(s).V,i?,(x(s))a;,(0))-(x,(0)p,(s).V,i?,(x(s)))] 
Jo 

= 0. (19) 
Thus, ( p6| ) reduces to 

{ix^{t) - X^iO)f) = 2K^t 

-(k,/3)2 f ds f ds' {pj{s - s')pfe(0)V,^j(x(s - s'))V.^fc(x(0))) . (20) 
Jo Jo 

This equation is however exact and can be evaluated to order by calculating the integral on 
its right hand side using the statistics for x and p in the absence of E, i.e. in the weak disorder 
limit which should become exact at high temperatures. Here p and x are purely Gaussian with 
correlation functions 



fe(^K(0))-|exp(^-^.)<5,fc, (21) 



and 



{{x,{t)-Xj{{))f)^25,jK,t. (22) 



This approximation yields 



{p,{s - s')pk{0)y^E,{^{s - s'))V,Sfc(x(0))) 

- ^ exp {-^{s - s')\ {V.E.i^is - s'))V.i?fc(x(0)))5,fe. (23) 



/9 V X 
Summing over the spatial indices then gives 



1^ ds' exp (^~^{s ~ s')^ (V,£;,(x(s-s'))V,£;,(x(0))).(24) 

The correlation function of the random field component in the above can be computed in terms of 
the correlation function A(x) of the random electrostatic field (p. 
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((v,i?,xGs-,s'))v.ii;,(x(o))) 



dk 



dk 



A{k)k'^ (exp(ik.x(s - s'))> 
^A(k)k4exp(-k2K,(s-s')) 



(25) 



where A(k) is the Fourier transform of A(x). Finally, we are left with a simple double integration 
over s and s', and from the large time behaviour we can extract Kg as 



kIi3x f dk A(k)k4 



(26) 



3. Analytic results in one dimension 



In the previous section, all the results we derived were for an arbitrary dimension D. To analyze our 
results to first order beyond the adiabatic limit, we will restrict our study to the one dimensional 
case where one can compute exactly the diffusion constant Kg for the effective Fokker-Planck 
equation (^2|). We apply the general results of ||l| to compute the effective diffusivity via a static 
problem: 



, , /dp ^ dV\ . 
\dx dx J 

where j is the current and k{x) is the spatially varying diffusivity in one dimension 

2' 



We can solve (E7h and find 



1 X 



'^ 2i dE{x) 
p 



dx 



Jo i^yy) 



dy 



(27) 



(28) 



(29) 



The diffusion constant Kg is then given as Kg \ = j- Where (.) is the average over the disorder 
in the random electrical field. This gives 



, dp\ Kej / . 

dx / L^co L \ 



and finally 



(30) 



(31) 



This result can also be obtained via a first passage time argument |23|. The electrical potential 
0(a;) is Gaussian and if we choose a correlator of the form A{x) — /(x^), where / is analytic at 
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X = 0, then we find that (^^^(O)) = -2/'(0), and {E'{0)E{0)) = 0. Therefore, E' and E are 
uncorrelated and we can write 

. - ^-'^)"' (32) 

In (corresponding to the adiabatic case here) it was shown that the dynamical transition could 
be identified, via an Arrhenius type argument, with the divergence of one or the other of the first 
two terms on the denominator of the equation above. We thus see that to first order beyond 
the adiabatic approximation the location of the transition temperature should be the same in the 
adiabatic limit and close to this limit. 

In the numerical simulations we will carry out we take the choice of correlation function f{u) — 
exp(— iy) for which wc find 

n^^n^-x'PY^'il-^). (33) 



4. Numerical simulations 



In this section we test our analytical predictions against numerical simulations of the Langevin 
equations (||) and (^) in one dimension. In our simulation, we set the diffusion constant without 
disorder to 1 for convenience and we used rescctled vcLria-bles pi — pi -yj f^p and Ei = Eiy/K^, so 
that both equations have the same time scale, 

^^/?.,^.g + y2,.. (34) 
5 = _ «p/3 f ^ - E,) + V2r,j,^ (35) 



dt - \x 

The coupled stochastic differential equations were integrated using second order Rungc-Kutta 
integration schemes developed in [|4| ^ and reviewed in [Q with a time step At — 0.001. In 
all simulations the effective diffusion constant for a given realization of the disorder was obtained 
by fitting the mean squared displacement averaged over 2000 particles at late times. The time of 
the simulation was chosen so that particles had typically diffused ten or so correlation lengths of 
the field. The fit of the average mean square displacement was done over the last half of the time 
of the simulation (to ensure that the mean squared displacement is well within the linear regime) 
by a non linear form: at + bt^ , with 6 < 1. The correction to the linear term has to be taken into 
account as it becomes more and more relevant close to the transition due to the slowing down of 
the dynamics [p^ . 

To generate the Gaussian field characterized by a correlator (^) , we used the technique presented 
in ||2^, |2^ . The process can be written in a general form 

K{x-x')r]{x')dx', (36) 
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where K is a kernel function and ri{x) is a Gaussian white noise. The linear form of (|3^), ensures 
that (f>{x) is a Gaussian process and the translationary invariance of the kernel K{x — x') ensures 
its stationarity. Now taking the Fourier transform defined by ^(w) — J^^^ (f>{x) exp{—iijjx),dx of 
we find 

4>{Lo)^K{oj)viLu). (37) 

If we now take the Fourier transform of the correlator (^ and the correlator {(j>{uj)(j){uj')) from (^7|), 
then by identification, we have A(w) = |i^((jj)p and thus from ( |3^ ) 

4>iLu) = ^/M^Mlu) (38) 

We can now create 4>{x) by sampling (I){uj) on a frequency mesh and taking the inverse fast Fourier 
transform (FFT). In jl^, |2j, the authors used a method due to Kraichnan which gives 

some finite size corrections for the diffusion constant close to the transition. However as we 
wish to precisely locate the transition and tests its dependence on Hp (and as we do not have 
analytical results for all parameter ranges) we use this FFT based method. Again we take a 
correlator of the form A(a::) = exp(— ^x^) and we average our sample over 2000 realizations of 
the field. In all our simulation we set x — ^- We measured the average value of Kg over the 
field with an error bar estimated from the standard deviation from sample to sample. The result 
in one dimension is compared in figure (^ for Kp = 10 with the analytical result ( ^3| ) and the 
Kubo formula (p6|). As expected, we are well in the adiabatic regime and ( |3^ ) is in a very good 
agreement with the numerical simulations. Moreover, we find a dynamical transition at /?c = 1 
where the diffusion constant becomes zero and below the critical temperature Tc, the diffusion 
turns out to be anomalous (see later). Finally, the agreement between (^3|) and the numerical 
results breaks down at high temperature, but fortunately in this regime the results agree with the 
weak disorder/high temperature analytical expression (p6|). We have thus a crossover between two 
analytical approximation which allows us to predict the behavior of Kg in the whole normal diffusion 
regime. 

Beyond the adiabatic approximation, that is to say when Kp is of the same order of Kx or smaller, 
we lack analytical results (except for weak disorder), thus wc must calculate Ke numerically. The 
results are plotted in figure (0) for Kp ~ 1 and figure (^) for Kp = 0.1 and compared with the weak 
disorder result of equation (P6|). In both cases, the transition still appears to occur at, or very 
close to, the same critical temperature Tc = I. Moreover, the diffusion constant clearly decreases 
monotonically with the value of Kp at fixed temperature. The low temperature phase, below T^, is 
characterized by an anomalous sub-diffusive behavior 

((x(^)-x(0))2)^c•^^ (39) 

where the exponent associated with the anomalous diffusion ly < 1/2. The subject of anomalous 
diffusion in disordered media has been extensively studied and a good review of it can be found in 
In the case of Kp ^ oo the authors in [|l^ were able to evaluate v by means of first passage 
calculation and replica trick and they found 
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Figure 1. Numerical value of Ke (stars) for Kp = 10, compared with analytical results ( [26 
(dashed line) and (pSl) (solid line). 




Figure 2. Numerical value of fte (plus) for Kp = 1, compared with analytical results (g6|) (solid 
line). 



We plot in figure ^ the value of the exponent v fitted from simulations with (^) for different values 



of Kp. The results agree relatively well with equation (|40|) for temperatures below but close to the 
transition temperature. However for very low temperature, the exponent is significantly larger than 
that predicted by equation (^0|). However we have verified that as the time of the simulation is 
increased that the measured exponent appears to decrease (indeed this was the case in ^| where 
analytical results and an effective trap model were available to identify the correct exponent). 
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Figure 3. Numerical value of Ke (plus) for Kp = 0.1, compared with analytical results (solid 
line). 



^0.3- 




Figure 4. Numerical value of u (cross) for Kp = 10, (stars) for Kp = 1 and (plus) for Kp = 0.1 
compared with analytical results (Uy) (solid line) in the low temperature phase. 



5. Conclusions and discussion 



We have studied the dynamics of a dipole diffusing in a random electrical field E derived from 
a quenched Gaussian potential. In the adiabatic limit (where the dipoles adapt very quickly to 
the local field), we showed that the coupled stochastic equation can be reduced, up to the order 
0{kx/kp), to an effective Langevin equation for a particle diffusing in an effective potential with a 
spatially varying and anisotropic local diffusivity Kij. In one dimension, we could compute exactly 
the diffusion constant of this process and we found a dynamical transition at finite temperature 
/3c = 1 with a crossover between a diffusive and a sub-diffusive regime. The validity of the effective 
Langevin equation ( |l^ breaks down for small /3, but a high temperature treatment allowed us to 
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compute a Kubo formula of Kg for a given value of Kp. We numerically checked our analytical 
predictions in one dimension finding good agreement and confirming the presence of the transition. 
We also ran simulations far from the adiabatic regime and we found strong indications that the 
transition remains at, or close to, the same (3c as for the adiabatic limit. Finally, we performed 
extensive numerical simulations in the low temperature phase to compute the anomalous exponent 
i/. We found that it does not depend on the relative values of diffusion constants of the particle's 
position and dipole moment and agrees rather well with a calculation made for the adiabatic case 




In [|2|, the authors showed that in higher dimensions, for the adiabatic case case Kp oo, that the 
same type of dynamical transition is also present. It would be interesting to investigate the dipole 
problem in higher dimensions to see if the conclusions of our current study remain valid. It would 
also be interesting to see if renormalization group type treatments or other approximation schemes 
could be developed in order to obtain analytical results in higher dimensions. 



Appendix A. Derivation of the effective long-time Langevin equation for x 



This appendix gives the detail of the calculation of the effective long-time Langevin equation for x. 
The first step is to formally integrate the equation (0) for the variable p to obtain 



where 



Piit) ^ Pi{0)exp{-^^t)+ [ dsexp ( 

X Jo \ X 



X Jo \ X 



(t-s) ) «p/3i?,(x(s))+j|-|G(i),(A.l) 



(A.2) 



is a colored noise. At late times the first term of equation (A.l) is exponentially suppressed and 
can thus be dropped giving 



/3k, 



p^(t)^ / dsexv\- — {t-s)\Kpl3E,{K{s)) + J—^Ut) 



X 



Substituting this into equation (0) yields 



2 X, 



dxi dE^ 
—r- = Pkx^^ 
dt oxj 



ds exp 



PKp ^ 

X 



Kp (3 



{t-s)] Kpl3E,{x{s)) + 



^X 

Kp (3 



(A.3) 



+ V2Ka;?7s,.(A. 



To procee d fu rther we will find an effective Fokker-Planck equation corresponding to the stochastic 
equation (A4). We follow the standard procedure of evaluating the change of an arbitrary function 
oix{t): 



, df{^{t)) 
' dt 



I 



d^Hp{^,t)f{^) 
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dxp(x, <)ffV(x) - (i/t/(x(t))> (A.5) 

where p(x, t) is the probabiUty density for the process x at time t. The operator H is the 
forward Fokker-Planck operator and ijt is its adjoint commonly called the backward Fokker-Planck 
operator. We now write using the Stratonovich prescription for white noise in stochastic calculus 
(where normal differentiation applies) to give 

= Ti+T2 + T3. (A.6) 

The first term is 

Ti = \/2^(VJ(x(0)?7x.), (A.7) 
which can be evaluated via Novikov's theorem as 

which is true for any Gaussian noise. Now for a white noise correlation, if we work with a 
symmetrised Dirac delta function corresponding to the Stratonovich prescription, the above integral 
picks up half the weight of the delta function. And we obtain 

The functional derivatives required in the above calculations can be evaluated using a path integral 
formalism as in pO| . However they can also be obtained directly from the stochastic equation as 
follows. Consider a general differential equation of the form 

dxi 
~dt 



W^ + B,,^,. (A. 10) 



In order to compute , I one can define the process 



_ &Xi{t) 

dr'^ 

-± = W. + B,,{^,+e,5{t-s)). (A.ll) 



Now from one of the standard definitions of the functional integral we have 
5x,{t) dxlit) 



6C,{s) de 



(A.12) 



Integration of equation ( A.ll ) in an infinitesimal interval about t = s yields 

x'{t+)^xit-) + B,,ej, (A.13) 

as long as the Wi do not depend on the derivative of x. We see that in the case considered here 
the Wi will be continuous at s — i and thus give no contribution to the functional derivative. This 
thus yields 

6xi{t) 



\s=t — J^ij 



B,j. (A. 14) 
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Using this result we obtain 
Sxj{t) 



which upon substitution in equation 

52 



) gives 



Tl — 



/(xO)- 



(A.15) 



(A.16) 



^ dxidxi 

This is the famihar Laplacian form arising for standard white noise. The second term T2 is given 

by 

I Kx /dEj df 



T2 



Once again using Novikov's theorem we find 



(A.17) 



T2 



^ y Kp dxk V dxi dxi 



' ds^-^{Qm{s)). 







(A.18) 



The correlation function of the Gaussian field C, is easily computed and is given for t>s by, 

(Ci WO(s)> = / du exp {-a{t + s- 2u)) , 



where a = Kp(3/x- At large t and s this becomes 

(OWO(s))-^cxp(-a(i-.)). 



(A.19) 



(A.20) 



We now consider the action of this correlation function as an operator on an arbitrary function g{s) 
via the following operator expansion 

rt 



ds exp {—as)) g{t — s) 



/W--4(t) + o(^) 

a at a 



and to leading order in 1/a we thus obtain 

■dE^ df\ 5xk{t), 



Kp dxk \dxi dxi J 5C,j{s) 



Now repeating the argument leading to equation (A.15) we find that 



5xk{t) 



dE^ 
dxk 



and thus 



T2 



,dEj d fdEj df 



' dxk dxk V dx. 



(A.21) 



(A.22) 



(A.23) 



(A.24) 



One thus sees that the noise C, is to leading order white noise interpreted via the Stratonovich 
prescription. The final term to evaluate is 



?3 = P^KxHpi 



dfi^{t))dEMt)) 



dxi 



dxi 



dsexp(— a(i — s))Ej{-x.{s))) 



(A.25) 
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and so one needs to compute terms of the form 

/ = (A(x(i)) / dsexp{-a{t - s))S(x(s))), 
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(A.26) 



bearing in mind that a is large. Now if p(x, t) is the probabiHty density for the position x and H 
the effective forward Fokker Planck operator, we may write the above as 



(ixS(x)p(x, s) exp(-(t - s)(a - H^))A{x). 



The density p obeys the Fokker Planck equation 



(A.27) 



(A.28) 



The integration over s can be evaluated by integrating by parts, we have 

/ ds p(x, s) exp(-(t - s)(a - H^)) = p(x, t) [l - exp(-i(a - H^))] {a - H^)-^ 
Jo 

ds exp(-(i - s)(a - H^)){a - H^)-\ (A.29) 

as 

Now because the eigenvalues of H must be negative or zero, we may neglect the second term in the 
square brackets on the right hand side above. Another integration by parts yields 

/ ds p(x, s) exp(-(i-s)(a-i7t)) = p(x, t){a-H^)-^-[Hp{x, t)] [a-H^Y^, (A.30) 

where we have used equation ( |A.28 ) in the second term of the right hand side and the operator H 
in this term only acts on p inside the square bracket. All other operators act on the right and we 
can now expand in powers of l/a to obtain 



[ ds S(x)p(x, s) exp(-(t - s)(a - H^))A(^) 
Jo 

p(x, t) 



A(X)B(X) , ^ , , , 

^ ' ^ ' + S(x)— A(x) - — ^(x)i? X) 



(A.31) 



Now putting all this together yields 

92 



■/(xO> 



— ) X 



,dEj d fdEj df 



This can be written in the form 



' dxk dxk V dxi dx 
-Ej -j- — -E^H 



OXi OXi KpP 



dxi dxi KpP dxi dx 



(i/t/(x)) = (ff„V(x)) 
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where 



2 '^x 

+ x — 



Hlf = nx[ 



dxk dxk V dxi dxi 



dxidxi 



dxi dxi •' 



This is the only term remaining strictly in the limit Kp 
particle in an effective potential V given by equation (|1 
we may write equation ( A.3S ) as 



dxi dxi 



dxi dxi 



(A.33) 



(A.34) 

oo and corresponds exactly to a Langevin 
. Now keeping the first term of 0{kx/kp) 



(i7t/(x)) = 



2 '^x 



dEj d fdEj df 



dxk dxk \ dxi dx. 



EM''"'' 



dxj dxj 



" dxj dxj 



(^o/(^)) '^2; ^Qxh V dxh dxi dxj 



^dE.dE, df dEi 

OXk OXi OXi OXk 



Finally the above may be written as 



(i?V(x)) = 



d 
dxk 



5/ 

dxi 



dV df . 
dxk dxi 



^). (A.35) 



(A.36) 



with Kij the spatially varying and anisotropic diffusivity tensor given by ( p!3[ ) and V the potential 
defined by (0). 
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